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随 着 信息 化 理论 对 军事 战争 的 深入 影响 ， 如 何 建 立 战场 环 。 的 数学 模型 实现 对 平台 的 聚 类 问题 ， 是 C2 组 织 决策 实体 配置 
过 程 中 主要 实现 的 内 容 呈 2 。 文 献 错误 ! 未 找到 引用 源 。 
(command and control, C2) 组 织 设计 中 的 关键 问题 二。 以 最 小 化 决策 实体 的 最 大 工作 负载 为 目标 函数 建立 了 数学 模 
通过 设计 作战 平台 与 决策 实体 之 间 的 指挥 控制 关系 ,使 。 型 ， 从 平台 、 任务 数量 的 角度 衡量 工作 负载 ; 文献 错误 ! 未 找到 


境 下 平台 资源 、 作 战 任务 、 决 策 实体 三 者 之 间 的 关系 ， 是 指挥 


控制 


引用 源 。 
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摘 要 : 为 解决 战场 C2 组 织 中 平台 资源 的 聚 类 问题 ， 首 先 定义 了 决策 实体 配置 过 程 中 各 个 要 素 ， 然 后 以 面向 作战 任 
务 的 形式 从 外 部 和 内 部 两 个 方面 定义 了 决策 实体 的 工作 负载 ， 以 工作 负载 的 均 方 根 (root mean square，RMS) 为 目标 
函数 建立 决策 实体 配置 模型 ,将 平台 调度 方案 作为 输入 信息 ,采用 基于 全 局 更 新 规则 的 蚁 群 算法 对 问题 模型 进行 求解 ， 
生成 符合 作战 需求 的 决策 实体 配置 方案 。 最 后 通过 仿真 算 例 对 所 提 方 法 进行 验证 和 分 析 ， 实 验 结 果 表 明 ， 该 方法 能 够 
得 到 较 好 的 平台 聚 类 方案 ， 在 聚 类 结果 和 收敛 速度 两 方面 上 与 传统 方法 相 比 均 具 有 优越 性 。 

关键 词 : 决策 实体 配置 ; 面向 作战 任务 ; 均 方 根 ; 蚁 群 算法 
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Solving method for decision-makers configuration problem based on 
global update rule of ant colony algorithm 
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Abstract: To solve the problem of clustering of platform resources in C2 organization, firstly, each element in process of 
decision-making entity configuration is defined, then the work load are defined by two aspects as inside workload and outside 
workload. The root mean square (RMS) of work load are established as the objective function of configuration model, and 
platform scheduling scheme is taken as the input information, global updating rule based ant colony algorithm is applied to 
solve this problem. Finally, the method are verified and analyzed by simulation example, results show that a greater clustering 
scheme are obtained by this method, and it has advantages over traditional methods in terms of clustering results and 
convergence speed. 
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日 所 错误 ! 未 找到 引用 源 。 错 误 ! 未 找到 引用 源 。 
及 o 


引言 ， 


对 


如 何 科 学 合理 的 对 工作 负载 进行 接 


述 和 测度 ， 并 建立 相应 


不 同 的 平台 在 相应 的 决策 实体 控制 下 能 更 好 更 快 地 完成 该 平台 ”引用 源 。 以 作战 任务 的 执行 时 间 对 工作 负载 进行 测度 ， 旭 
所 分 配 到 的 作战 任务 ， 提 高 任务 执行 效果 和 执行 质量 ， 就 是 平 ”以 工作 负载 的 最 小 RMS 值 为 目标 函数 的 数学 模型 ， 文 献 错误 ! 
台 资 源 的 聚 类 问题 ， 也 称 为 C2 组 织 设计 中 的 决策 实体 配置 问 ” 未 找到 引用 源 。 中 定义 了 指控 复杂 度 的 概念 ， 从 不 同 平台 对 作 
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录用 稿 陈 冠 宇 ， 等 : 基于 全 局 更 新 规划 皮革 弦 和 有 决策 实体 配置 问题 求解 方法 
战 任务 的 协作 情况 建立 最 小 化 工作 负载 的 问题 模型 ; 文献 错误 ! e) 在 执行 作战 任务 中 ， 定 义 决策 实 体 集 合 


未 找到 引用 源 。 针 对 传统 层次 聚 类 方法 采用 的 贪心 策略 容易 陷 
入 局 部 最 优 的 现象 ， 提 出 备 选 解 策略 用 于 选 出 更 加 满足 整体 需 
求 的 最 优 解 ; 针对 多 目标 的 优化 问题 , 文献 错误 ! 未 找到 引用 源 。 ” 义 所 有 的 作战 平台 集合 S$” = { 有 ,局 …Pe}，NP 表示 作战 平台 


i 


DM ={DM1,DM,,…, DMwp} ， 其 中 ，ND 为 决策 实体 数量 ， 定 


提出 对 多 目标 问题 模型 生成 一 系列 前 沿 解 ， 并 通过 一 定 的 筛选 ” 数量 ;定义 所 有 的 作战 任务 集合 5S” ={T,,…,Tyr} ，NT 表示 
方法 得 到 符合 要 求 的 平台 聚 类 结果 。 所 有 任务 数量 。 

本 文 在 以 上 文献 研究 内 容 的 基础 上 ， 进 一 步 讨 论 针 对 作战 ” 1.2 约束 条 件 分 析 
任务 外 部 工作 负载 和 内 部 工作 负载 结合 描述 决策 实体 工作 负载 兵力 组 织 的 结构 可 分 为 决策 层 、 任 务 层 和 平台 层 ， 这 
测度 方法 , 建立 以 决策 实体 工作 负载 RMS 值 为 目标 函数 的 决策 。 中 。 
实体 配置 问题 模型 ， 利 用 基于 全 局 更 新 规则 的 蚁 群 算法 对 问题 a) 决策 实体 和 作战 平台 是 一 个 一 对 多 的 关系 。 一 个 决策 实 
模型 进行 求解 ， 最 后 通过 仿真 算 例 对 算法 最 佳 参 数 设 定 进行 讨 ” 体 可 以 指挥 控制 多 个 作战 平台 ， 而 一 个 作战 平台 只 能 隶属 于 一 


论 ， 并 对 算法 的 适用 性 与 优越 性 进行 验证 与 分 析 ， 实 验 结果 表 ”个 决策 实体 。 即 
明 选 择 合适 的 算法 参数 能 够 得 到 科学 合理 的 决策 实体 配置 方 Pm 1 ( j=12,., NP) 
案 ， 并 相 较 于 传统 蚁 群 算法 、 自 适应 蚁 群 算法 ， 该 算法 具有 较 nl 
快 的 收敛 速度 且 能 得 到 较 好 的 聚 类 结果 。 b) 由 变量 描述 可 知 , 当 且 仅 当 xuw =1 且 入 =1 时 ，yw =1l。 
题 措 述 及 建 模 因此 存在 关系 
Vim S No (5) 
1.1 变量 定义 nd es ND,mzn) 
战场 环境 中 的 所 有 实体 主要 被 分 为 三 类 ， 分 别 为 决策 实体 同 理 可 得 
(CDM )、 平 台 ( Pf; ) 以 及 作战 任务 〈 五 )。 平 台 聚 类 问题 ， mi S Xo (0) 
即 在 任务 与 平台 之 间 的 分 配 关系 产生 的 前 提 下 ， 实 现 平 台 的 聚 (m=1,2,..., ND;n=1,2,..., ND,mzn) 
类 并 将 其 分 配给 不 同 决 策 实体 进行 指挥 控制 的 过 程 。 针 对 本 文 为 为 大 三 四 xm ， 所 以 
问题 ， 首 先 对 相关 变量 作 如 下 定义 : Vim < Oy am (7) 
a) 平台 与 作战 任务 的 匹配 关系 用 矩阵 (mn=1,2,..., ND,m 关 7i= 上 2…,NT;7 =12…,NP) 同 理 可 
N=( 四 )，，G=12.…NTJ=L2.…NP) 表 示 。 若 几 =1 ， 则 表 人 人 加 
示 任 务 卫 需要 平台 已 执行 ， 若 由 =0 ， 则 表示 任务 五 不 需要 平 c) 每 一 个 决策 实体 必须 参与 指挥 控制 平台 执行 该 平台 相应 
台 咏 执行 。 该 矩阵 为 模型 输入 已 知 信息 。 分 配 到 的 作战 任务 ， 即 每 个 决策 实体 指挥 控制 的 平台 数 必须 至 
b ) 平台 与 决策 实体 的 匹配 关系 用 矩阵 少 为 1 个 : 
M=(m) =1L2…ND;7 =12…NP) 表示 。 若 wm =1 则 表 pi (n=1,2,..., ND) (9) 


示 决 策 实体 DM 指挥 控制 平台 了 ， 若 m=0 则 表示 决策 实体 。 ”1.3 问题 模型 建立 


DM 不 对 平台 PD 进行 指挥 控制 。 此 和 矩 阵 即 为 模型 最 终 输出 信 ”1.3.1 工作 负载 的 定义 
息 。 工作 负载 表明 了 每 个 决策 实体 所 指挥 控制 的 平台 执行 任务 
c) 任务 与 决策 实体 的 关系 用 和 矩阵 的 情况 。 考 虑 到 作战 任务 的 完成 时 间 有 差异 ， 并 且 对 不 同 的 决 
i= =X) vpxwr (Nn =D >. ND;i=1,2,.…,NT) 表示 : 策 实体 带 来 的 工作 负载 量 不 同 ， 本 文选 择 以 任务 执行 时 间作 为 
1， 表 示 决 策 实体 DMn 执 行 任务 Ti 工作 负载 测度 并 从 内 部 负载 和 外 部 负载 两 个 方面 
.二 表示 决策 实体 DMn 不 执行 任务 Ti - 进行 描述 。 
i 可 由 和 mi 的 关系 推导 得 出 ， 即 车 对 于 平台 PD， 其 1) 内 部 负载 1(m) 定义 及 测度 方法 
中 =1 且 m=1， 则 。 决策 实体 内 部 负载 1(n) 表示 决策 实体 指挥 控制 隶属 于 自 
d) 决策 实体 之 间 的 协作 关系 用 向 量 y, 表示 : 己 的 作战 平台 完成 作战 任务 所 产生 的 负载 。7(n) 数值 上 等 于 决 


策 实体 控制 的 所 有 平台 的 执行 作战 任务 的 时 间 之 和 : 


0 表示 DAM ,和 DAM ,在 任务 T 上 存在 协作 
ke DM, DM T 上 不 存在 协 Np 
0， 表 示 DM,, 和 DM ,在 任务 上 不 存在 协作 100) -Pm, .sn 而 
因此 ， 以 上 描述 可 知 : j=1 
Vom = Xo Xs G) ”其 中 : 5(7) 表示 平台 执行 任务 的 工作 负载 , 定义 为 作战 平台 
(mn =1,2,.., NDym #1:i =1,2,.., NT 执行 任务 的 时 间 之 和 |: 


(2) 


可 


作 负 载 的 不 同 将 对 配 
的 平台 调度 方案 ， 无 论 聚 类 结果 如 何 变化 ， 所 有 决策 实体 的 内 
部 负载 之 和 是 保持 不 变 的 ， 但 不 同 的 聚 类 方案 决策 实体 的 外 部 
负载 之 和 却 不 同 。 因 此 可 以 通过 选择 合适 的 聚 类 方案 使 得 所 有 
决策 实体 的 外 部 负载 之 和 最 小 。 本 文选 择 用 每 一 个 决策 实体 平 
均 工 作 负载 最 小 、 同 时 决策 实体 之 间 负 载 差异 最 小 作为 聚 类 模 
型 的 目标 函数 ， 其 中 ， 决 策 实体 的 平均 工作 负载 表示 为 
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Ina 


ch 
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5())=20,.t GD) 
1 


其 中 ， 厂 表示 任务 五 的 执行 时 间 。 


因此 ， 综 上 式 (10) (11) 两 式 可 知 : 


[03 CD 5 (12) 


2) 外 部 负载 £(n) 定义 及 测度 方法 
决策 实体 外 部 负载 表示 两 个 不 同 决策 实体 协作 完成 作战 任 


务 而 产生 的 工作 负载 。E(n) 数值 上 等 于 决策 实体 PM, 与 其 他 
所 有 决策 实体 控制 的 平台 协作 执行 任务 的 工作 负载 之 和 。 类 比 
决策 实体 内 部 负载 的 定义 方法 将 E(n) 定义 为 


ND 


EW= >, Dy (13) 
L 


m=1,mz#n i= 


3) 决策 实体 DM 总 工作 负载 CW(n) 测度 方法 
综合 式 (12) (13)， 决 策 实体 总 工作 负载 CW(n) 的 定义 为 

CW(n)=a:1(n)+b:E(n) (14) 
其 中 : 权重 系数 a 和 4 分别 为 内 部 负载 和 外 部 负载 对 总 工 


作 负 载 的 影响 ，a+b=1。 
1.3.2 目标 函数 构造 


于 工作 负载 从 内 部 负载 和 外 部 负载 两 个 方面 进行 描述 
的 聚 类 方案 中 决策 实体 指挥 和 控制 作战 平台 完成 任务 的 工 
方案 产生 不 同 的 结果 。 对 于 某 一 个 给 定 


1 D 
/= 二 2 ,CW(m) (15) 


m=1 


决策 实体 之 间 的 差异 由 方差 表示 ， 即 


=ICWom -4 (10) 


m=1 


因此 ， 用 均 方 根 RMS 值 对 决策 实体 的 工作 负载 进行 综合 


者 述 ， 即 


| 1 © | 1 2 
RMS=| 2 CW = EW +2 一 2 人 
m=1 m=1 


D 
+ > (Cwony + 47 —24: CW(m)) 


m=1 


= to (17) 


陈 冠 宇 ， 等 : 基于 全 局 更 新 规则 蚁 群 优化 的 


1.3.3 数学 模型 建立 
式 (14) 可 知 ， 决 策 实体 的 工作 负载 将 通过 均值 和 方差 两 
因此 ， 在 对 问题 模型 的 约束 条 件 分 析 及 目标 函 


方面 进行 描述 ， 
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的 决策 实体 配置 问题 求解 方法 


数 构造 的 基础 上 ， 本 文 建立 的 决策 实体 配置 问题 数学 模型 表示 


为 


min RMS= SY cwon (18) 
m=1 


ND 
st. Dum, =1, j=1,2,...,NP 


n=1 


NP 
Dm >1 i=1,2,..., NT; 
n=1 


yumi < Oj m 


m=1,2,..., ND; n=1,2,...,ND 


Yi < Oy “Mj, m=1,2,..., ND; n=1,2,...,ND 


2 ”基于 全 局 更 新 规则 蚁 群 算法 的 决策 实体 配置 问题 


求解 


决策 实体 配置 草 


是 实现 平台 聚 类 的 过 程 ， 即 把 不 同 的 平台 


分 配给 一 定数 量 的 六 
定 的 作战 任务 。 在 3 
同时 得 到 符合 问题 模型 
体现 了 所 设计 的 


策 实 体 ， 由 决策 实体 控制 各 个 平台 执行 既 
F 台 聚 类 过 程 中 ， 要 在 满足 问题 约束 条 件 的 


标 函 数 的 最 优 解 ， 因 此 平台 聚 类 方案 
F 台 与 决策 实体 之 间 配 置 关系 的 合理 性 。 


决策 实体 的 配置 问题 本 质 上 是 组 合 优化 问题 ， 可 以 采用 智 


能 算法 对 问题 进 


ACO), 作为 一 种 


法 ， 已 经 被 广泛 


9 、 流 水 车 间 调 
选 址 第 误 ! 未 的 到 引用 源 。 和 无 


着 较 好 的 表现 ， 


该 算法 按照 
j 蚂 蚁 在 所 通过 


Sig 


迭代 逐步 收敛 到 所 求 和 
以 简要 描述 为 : 月 
同 蚂蚁 行走 的 路 径 对 应 不 同 


行 求 解 。 蚁 群 算法 (ant colony optimization， 
具有 搜索 能 力 强 、 收敛 速度 快 等 特点 的 智能 
于 求解 组 合 优化 问题 。 在 路 径 规 划 和 下 和 9 
度 学 各 林 到 引用 浙 ,、 最 优 解 搜寻 千克 到 引用 滨 ， 建筑 
人 机 作战 目标 分 配 下 有, 等 问题 上 有 
此 蚁 群 算法 可 以 作为 解决 平台 聚 类 过 程 中 寻 
求 最 优 聚 类 方案 的 一 种 方法 。 
思想 ,使 用 正 反 馈 和 信息 素 挥 发 机 制 ， 利 
上 释放 的 信息 素 的 诱导 作用 ,通过 多 次 
题 的 全 局 最 优 解 。 算 法 的 核心 思想 可 
会 在 经 过 的 路 径 上 留 下 信息 素 ， 不 


的 蚂蚁 会 更 早 的 到 达 


的 信息 素 轨迹 。 由 于 在 最 短路 径 上 


的 地 ， 从 而 增加 其 在 该 路 径 上 来 往 的 频 


率 ， 因 此 该 路 径 上 信息 素 密度 更 大 。 其 他 蚂蚁 在 寻 径 的 过 程 中 


会 选择 信息 素 强 的 路 径 进 


行 移动 ， 因 此 ， 此 路 径 上 的 信息 素 强 


度 会 不 断 县 加 。 通 过 这 种 了 


到 食物 源 的 最 优 路 径 。 
下 面 将 针对 决策 实体 的 配置 问题 , 结合 文献 错误 ! 未 找到 引 


E 反 馈 机 制 使 蚂蚁 最 终 能 够 找到 蚁 六 


用 源 。， 对 蚁 群 算法 在 实现 平台 聚 类 过 程 中 算法 各 个 组 成 要 素 


进行 相应 的 描述 
a) 构建 图 。 


的 随机 构建 过 程 。 划 
包含 所 有 平台 和 决策 实体 。 


蚂蚁 将 在 构建 矩阵 C 上 随机 游 走 ， 来 建立 解 


工 完全 连接 集合 C 中 所 元素 ， C 中 


个 可 行 解 就 表示 一 个 由 ND 个 决 


策 实 体 和 NP 个 平台 以 (mm 了 的 配对 形式 建立 的 矩阵 ， 若 (mm 放 
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等 于 1， 则 表示 平台 已 分 配给 决策 实体 DM,, 进行 指挥 控制 。 
分 配 过 程 中 ， 一 个 平台 只 能 分 配给 一 个 决策 实体 ， 而 一 个 决策 
实体 可 以 指挥 控 台 。 

b) 信息 素 。 信 息 素 5 表明 平台 P 到 决策 实体 DM 的 其 
望 度 ， 它 编码 了 有 关 蚂 蚁 进行 路 径 搜索 的 过 程 ， 由 蚂蚁 对 其 进 
行 更 新 。 

c) 信息 素 全 局 更 新 。 为 了 避免 元 余 的 信息 素 造成 过 量 的 信 
息 堆积 ， 每 只 蚂蚁 完成 对 NP 个 平台 的 聚 类 后 将 淡化 旧 信息 ， 
强化 新 信息 。 因 此 ， 在 t+n 时 刻 ， 平 台 与 决策 实体 之 间 的 配对 
(mj 将 按 全 局 更 新 规则 进行 更 新 : 

T(t+n)=0 pT, (t+Ar, (tt+n) (19) 


趾 多 个 平台 


其 中 : P 为 挥发 系数 ， 表 示 信 息 素 随 时 间 推 移 的 衰减 程度 ; 
Azw 表示 信息 素 增 量 。 
二 (20) 
k=1 


Az(4bt+n) 表示 第 只 蚂蚁 在 (4t+n) 时刻 留 在 路 径 
(m,j) 上 的 信息 素 量 ， 具 体 数 值 根据 蚂蚁 的 表现 而 定 。 计 算 公 
式 为 


0 了 HA 、 人 中 
» 家 大 4 曙 经 十 全 站 
8 若 第 kK 只 蚂蚁 经 过 咯 径 (m 让 人 


0, 其 他 
其 中 : Q 表示 控制 常数 ，RMSi 表示 每 只 
方案 的 工作 负载 RMS 值 : 


ND 
RMS = 2 CW(n)’ = ye +0o” (22) 


d) 结束 规划。 如果 所 有 平台 都 已 分 配 完 
次 数 超过 某 个 迭代 次 数 限制 ， 则 结束 。 

综 上 所 述 ， 基 于 全 局 更 新 规则 蚁 群 算法 的 决策 实体 配置 过 
程 步 又 如 下 : 

a) 初 始 化 所 有 参数 。 设 置 蚂蚁 数量 为 N ， 最 大 循环 次 数 为 
Nsop。， 初 始 化 每 个 配对 (m7 的 信息 素 Tw(0) ， 令 初始 时 刻 
ATm(0)=0 。W 只 蚂蚁 随机 产生 决策 实体 与 平台 之 间 的 配置 方 
案 C ， 对 方案 的 可 行 性 进行 验证 ， 剔 除 不 可 行 的 配置 方案 ; 

b) 计 算 不 同 的 决策 实体 配置 方案 中 决策 实体 的 工作 负载 
RMS 值 ; 


Th (t,t+n)= 


Ty 


蚂蚁 表示 的 分 配 


， 或 者 算法 迭 代 
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昌 计 算 信息 素 增 量 人 A 态 (41+m) ,求解 信息 素 rv(r+ 癌 ， 


根据 信息 素 计算 下 式 的 转移 概率 
pr 2 Twn 
” > tn) (23) 


seallowed,. 


dd) 按照 转移 概率 进行 聚 类 方案 选择 ， 


记录 当前 最 优 分 配方 


e) 判 断 循 环 次 数 是 否 小 于 最 大 循环 次 数 ， 若 小 于 则 保存 最 
优 个 体 ， 返 回执 行 步 又 2， 若 大 于 则 输出 最 后 结果 。 
基于 蚁 群 算法 的 平台 聚 类 过 程 如 图 1 所 示 。 
(开始 ) 


Y 
初始 化 算法 参数 


输入 平台 -任务 分 配 关系 


(um 


基于 全 局 更 新 规则 的 下 好 
算法 产生 聚 类 方案 


计算 转移 概率 ， 根 据 概率 


一 二 
ET 产生 新 的 聚 类 方案 
立 3 
一 是 >、 更 新 信息 素 


NC<NC_MAX? 


输出 决策 实体 配置 方案 
六 


图 1 基于 全 局 更 新 规则 蚁 群 算法 的 聚 类 过 程 


3 ” 算 例 分 析 


本 文 以 一 个 联合 作战 作战 想 定 为 仿真 算 例 ， 对 本 文 所 提 方 
法 进行 仿真 验证 分 析 。 在 该 案例 中 ， 共 有 18 个 任务 、20 个 平 
台 ， 其 中 平台 资源 能 力 与 任务 资源 需求 的 种 类 为 8 种 ， 任 务 属 
性 数据 和 平台 属性 数据 参见 文献 错误 ! 未 找到 引用 源 。。 其 中 ， 
任务 -平台 间 的 分 配 关系 用 平台 调度 方案 甘 特 图 进行 表示 , 该 分 
配 关系 为 决策 实体 配置 问题 的 已 知 输入 信息 ， 如 图 2 所 示 。 


201806.00120v1 


chinaXiv 


http://www.arocmag.com/article/02-2019-10-014.html 


仿真 实验 1 在 基于 全 局 更 新 规则 的 蚁 群 算法 中 
涉及 到 的 参数 较 多 ， 且 不 同 参数 设 


类 结果 产生 


ChinaXiv 命 作 期 天 


算 机 应 用 研究 


于 算法 


较 大 影响 ， 因 此 


的 影响 进行 


比较 分 析 ， 确 定 


的 结果 将 对 最 终 的 平台 聚 
首先 针对 不 同 算法 参数 对 聚 类 结果 


解决 该 问题 的 最 佳 参 数 设置 问题 。 


1) 决策 实体 数量 ND 对 聚 类 结果 的 影响 


算法 初 


始 化 参数 设 


NP=20、 任 务 个 数 NT =18 


素 残 


系数 Rio=0.8 、 信 


如 


下 : 蚂蚁 数量 M =25 、 平 台 个 数 
、 最 大 循环 次 数 NC =50 ， 信 息 


a=0.5 、 外 部 负载 权重 b= 
望 值 启发 式 因 子 6=3 。 令 决策 实体 数量 分 别 为 ND=2,3,…,10， 
工作 负载 RMS 值 随 ND 变化 情况 如 图 3 所 示 。 


最 小 RMS 值 


500 


400 


300 


200 


2 4 


| 


息 素 总 


决策 实体 个 数 ND 

图 3 RMS 值 随 ND 个 数 变化 
3 可 知 ,决策 实体 数量 将 影响 目标 函数 值 的 变化 勤快 ， 
一 般 决 策 实体 的 个 数 越 少 ， 


EQ=100 、 内 部 负载 权 导 


图 2 平台 调度 方案 甘 特 图 


时 间 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 
Pl 7T17 TI12 

Pp2 Tl T2 

Ls | 下 

Pp4 T7 TI1 TIS T16 
ps ， TI7 T8 T15 

Pp6 Tl 

P7 4 TI1 T14 

pg Tl T2 

Pp9 13 17 T8 T12 

P10 7T6 17 T9 T13 TIS TI10 T16 
PIL TS T13 

Pl2 Ts m9 | T13 | 

Pl3 Ts Tg TI10 

pl4 5 T18 TI2 

Pl5 TI7 T2 T14 

P16 Tl T2 T14 

P17 T18 T13 

P18 16 T18 T16 
pl9 6 T7 TIS 

P20 TI17 | TS T16 


造成 决策 实体 成 本 增加 。 线 变化 趋势 可 以 看 出 ， 当 决策 实 
体 在 4-6 个 之 间 时 ， 曲 线 的 斜率 变化 最 快 ， 因 此 决策 实体 个 数 
应 该 选择 4-6 个 。 综 合 考虑 ， 选 择 决 策 实体 数量 ND =5 作为 后 
续 聚 类 过 程 的 参数 值 。 

2) 信息 素 残留 系数 CD 对 聚 类 结果 的 影响 

言 息 素 挥发 度 是 指 蚂蚁 在 走 过 的 路 径 上 留 下 的 信息 素 的 挥 
发 程度 。 在 算法 中 以 1-p 为 信息 素 挥发 系数 ,， 即 以 P 为 信息 素 


0.5 、 信 息 启 发 式 因 子 ga=0.8、 期 


1 线 图 


目标 函数 RMS 的 值 就 越 大 , 表明 当 


收 稿 日 期 : 2018-03-26; 修 回 日 期 : 2018-05-10 


决策 实体 数量 少时 ， 每 个 决策 实体 的 负载 越 大 。 但 是 并 不 是 决 
策 实体 的 个 数 越 多 ， 聚 类 方案 就 越 好 。 受 到 作战 实际 
决策 实体 个 数 是 有 限 的 ， 在 一 次 战役 任务 中 ， 过 多 决策 实体 将 


的 约束 ， 
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计算 机 可 视 化 . 


EE 残留 系数 。 信 息 素 挥发 系数 直接 影响 到 蚁 群 算法 的 收 和 敛 速度 和 
全 局 搜索 能 力 。 由 于 信息 挥发 系数 的 存在 ， 蚂 蚁 没有 走 过 的 路 


径 上 的 信息 素 会 逐渐 挥发 直至 减少 到 0。 当 


言 恩 素 挥 发 系数 
1-p 过 大 时 , 信息 素 挥 发 过 快 , 一 些 蚂蚁 未 经 过 路 径 不 会 得 到 
选择 ,降低 了 算法 的 全 局 搜索 能 力 ; 而 当 信 息 素 挥 发 系数 1-p 
过 小 时 ， 虽 然 能 提高 算法 的 全 局 搜索 能 力 ， 但 是 算法 收敛 速度 
会 降低 。 信 息 素 残 留 系数 对 算法 性 能 的 影响 如 表 1 所 示 


表 1 信息 素 残 留 系数 对 算法 性 能 的 影响 
残留 系数 PP 最 优 RMS 值 收敛 代数 
0.1 288.6357 20 
0.3 288.7796 15 
0.5 288.2410 11 
0.7 275.6084 15 
0.8 276.3512 16 
0.9 267.0019 21 
1 270.2591 19 
言 息 素 残 留 系数 P 与 算法 输出 RMS 的 之 间 关 系 如 图 4 所 


不 。 
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录用 稿 陈 冠 字 ， 等 : 基于 实体 配置 问题 求解 方法 
300 后 续 聚 类 过 程 的 参数 值 。 
4) 启发 因子 w 的 选择 
信息 启发 式 因 子 & 反映 蚂蚁 在 运动 过 程 中 所 积累 的 信息 
量 在 指导 蚁 群 搜索 中 的 相对 重要 程度 ,期 望 值 启发 式 因 子 B 反 
起 240 映 昌 蚁 在 运动 过 程 中 启发 信息 在 指导 蚁 群 搜索 中 的 相对 重要 程 
220 度 。 期望 值 启发 式 因 子 B 的 值 越 大 ,蚂蚁 在 某 个 局 部 点 上 选择 局 
zooL L L L | 部 最 短路 径 的 可 能 性 越 大 ,算法 收敛 速度 得 以 加 快 ,但 蚁 群 在 最 
信息 素 残留 系数 p 优 路 径 的 搜索 过 程 中 随机 性 减弱 ,易于 陷入 局 部 最 优 ; 而 信息 局 
图 4 信息 素 挥发 系数 与 输出 解 间 的 关系 发 式 因子 & 的 大 小 则 反映 了 蚁 群 在 路 径 搜索 中 随机 性 因素 作 
信息 素 残 留 系数 P 与 算法 收敛 迭代 次 数 之 间 的 关系 如 图 5 用 的 强度 ,其 值 越 大 ,蚂蚁 选择 以 前 走 过 的 路 径 的 可 能 性 越 大 , 搜 
所 示 索 的 随机 性 减弱 , 当 值 过 大 也 会 使 蚂蚁 的 搜索 过 早 陷于 局 部 最 
天 优 。 启发 因子 & 和 对 输出 结果 及 算法 收敛 性 的 影响 如 表 3 所 
20| 不 。 
表 3 启发 因子 & 和 对 输出 结果 及 算法 收敛 性 的 影响 
和 六 息 启发 因子 a 期 望 值 启发 因子 。 最 优 RMS 值 ”算法 收敛 代数 
EF 加 0.1 0.1 281.5496 80 
0.3 0.5 282.9456 40 
0.5 1.5 266.1665 25 
W 0.8 3 265.3575 18 
图 5 信息 素 残留 系数 p 与 算法 收敛 迭代 次 数 之 间 的 关系 1 5 266.4062 12 
线 趋势 可 以 看 出 ,算法 在 p=0.7 一 0.9 间 时 ,输出 RMS 3 8 272.4683 12 
最 小 。 因 此 本 设计 综合 其 收敛 速度 和 全 局 搜索 能 力 等 因素 。 综 5 10 296.6914 6 
合 考 虑 ， 选 择 信息 素 挥 发 系数 1-p=0.3 ， 即 信息 素 残 余 系数 由 表 3 可 知 当 信息 启发 因子 w<=0.3-3， 同 时 期 望 值 启发 
D=0.7 作 为 后 续 聚 类 过 程 的 参数 值 。 因子 6=3 左 右 时 算法 的 输出 结果 较为 优越 ， 且 算法 收敛 速度 
3) 蚂蚁 数量 M 对 聚 类 结果 的 影响 较 快 ， 因 此 设置 参数 w =0.8、Z =3 作 为 后 续 聚 类 过 程 的 参数 
蚁 群 算法 中 ， 蚂 蚁 的 数量 的 多 少 直接 影响 着 群体 表现 出 的 值 。 
能 力 。 对 于 单个 蚂蚁 而 言 完成 一 次 分 配 即 为 一 个 有 效 解 ， 所 有 仿真 实验 2 利用 实验 1 中 的 结果 ， 对 算法 中 的 参数 进行 
蚂蚁 完成 一 次 分 配 可 构成 一 个 解 集 。 显 然 ， 蚂 蚁 的 数量 越 多 ， 初始 化 设置 : 蚂蚁 数量 M =25 、 决 策 实体 数量 ND =5、 平 台 
解 集 越 大 ， 算 法 的 全 局 搜索 能 力 就 越 强 ， 但 是 算法 的 收敛 速度 。 个 数 NP=20、 任 务 个 数 NT =18 、 最 大 循环 次 数 NC =50 、 
就 越 慢 。 反 之 ， 曲 蚁 的 数量 越 少 ， 解 集 越 小 ， 算 法 的 全 局 搜索 言 息 素 残留 系数 Rho =0.8 、 信 息 素 总 量 Q=100 、 内 部 负载 权 
能 力 就 越 弱 , 但 是 算法 的 收敛 速度 就 越 快 。 当 昌 蚁 数量 过 小 时 ， 重 a=0.5 、 外 部 负载 权重 5=0.5、 信 息 启 发 式 因 子 & =0.8、 
虽然 提高 了 算法 的 收敛 速度 ,但 是 却 容易 使 算法 陷入 局 部 最 优 ， 期 望 值 启发 式 因子 =3 。 利 用 基于 全 局 更 新 规则 的 蚁 群 算法 ， 
全 局 搜索 效能 降低 ， 容 易 使 算法 过 早 停滞 现象 。 蚂 蚁 数量 M 对 生成 以 图 1 中 所 示 的 平台 -任务 分 配 关系 为 输入 信息 的 平台 聚 
对 聚 类 结果 及 算法 收敛 情况 的 影响 如 表 2 所 示 。 类 方案 ， 其 聚 类 结果 如 表 4 所 示 。 
表 2 蚂蚁 数量 M 与 RMS 值 、 算 法 迭代 次 数 关 系 表 表 4 基于 全 局 更 新 规则 蚁 群 算法 的 平台 聚 类 结果 
蚂蚁 数量 M 最 优 RMS 值 算法 收敛 代数 决策 实体 平台 DM 工作 负载 ”工作 负载 RMS 值 RMS 方差 
20 269.6683 18 DM' Pi, Pa, Ps, Ps, Po 260 
25 256.4048 20 DhM P,Ps, Ps, Ps Pis.Pis 350 
30 278.9426 22 DM Pr, PH 130 256.2421 80.1873 
35 273.5055 25 DMa Ps, Pu, Pu, Pus, Ps 220 
40 278.3597 30 DMs Pio Pis, Ps 270 
50 259.5824 34 算法 收敛 曲线 如 图 6 所 示 。 
80 272.2179 40 
表 2 可 知 ， 当 蚂蚁 数量 过 多 时 ， 收 敛 速度 降低 的 同时 会 
影响 算法 输出 的 聚 类 结果 。 由 图 可 以 看 出 选择 蚂蚁 只 数 M = 25 
时 ,输出 结果 较 小 、 迭 代 次 数 较 少 ， 因 此 设置 参数 M = 25 作为 
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360 采用 全 局 更 新 规则 的 蚁 群 算法 。 进 一 步 证 实 了 蚊 群 算法 在 解决 
340 平台 聚 类 问题 上 的 可 行 性 和 优越 性 。 
320 对 该 问题 分 别 采 用 全 局 更 新 规则 的 蚁 群 算法 、 自 适应 蚁 群 
00 算法 、 普 通 蚁 群 算法 和 量子 遗传 算法 (quantum genetic algorithm， 
宣 ，。 QGA) 进 行 求解 , 对 算法 收敛 情况 进行 进一步 比较 分 析 , 各 算法 
加 收敛 曲线 如 图 7 所 示 。 
渍 - _ 算法 对 比 图 _ _ 
由 算法 法 代 次 数 NC 区 - % 趾 一 -一 采用 全 周章 新 规则 的 数 群 法 | 
a 普通 蚁 群 算法 
图 6 蚊 群 算法 收敛 曲线 与 迭代 次 数 E ol 生子 法 估算 法 
图 6 及 表 4 可 知 ， 基 于 全 局 更 新 规则 的 蚁 群 算法 在 解决 ao 
平台 聚 类 问题 时 具有 可 行 性 。 从 算法 输出 结果 可 以 看 出 ， 决 策 a0 
实体 的 工作 负载 平均 值 为 246， 工 作 负载 RMS 值 为 256.2421， | 
决策 实体 工作 负载 方差 为 80.1873。 从 算法 收敛 曲线 可 以 看 出 ， 260|- 
算法 在 迭代 至 33 代 的 时 候 收敛 。 因此 , 蚁 群 算法 在 解决 解决 平 ~ 
台 聚 类 问题 上 具有 较 好 的 性 能 ， 验 证 了 该 方法 的 可 行 性 和 适 | 0 
性 。 图 7 不 同 算法 收敛 代数 比较 图 
仿真 实验 3 为 了 进一步 验证 基于 全 局 更 新 规则 的 蚁 群 算法 由 图 7 中 各 算法 收敛 曲线 趋势 可 以 直观 看 出 ， 采 用 全 局 更 


的 优越 性 ， 首 先 将 本 文 所 提 算 法 与 层次 聚 类 算法 、 自 适应 蚁 群 。 新 策略 的 蚁 群 算法 收敛 效果 与 采用 自 适应 蚁 群 算法 的 收敛 效果 
算法 进行 对 比分 析 。 算 法 中 的 关键 参数 设置 为 : 蚂蚁 数量 类 似 ， 全 局 更 新 策略 的 蚁 群 算法 的 收敛 效果 和 达到 收敛 的 迭代 
M =25 、 决 策 实体 数量 ND =5、 平 台 个 数 NP=20、 任 务 个 数 ” 数 均 优 于 量子 遗传 算法 ， 因 此 ， 进 一 步 验证 了 采用 全 局 更 新 策 
NT=18 、 最 大 循环 次 数 NC =50 、 信 息 素 残留 系数 略 的 蚁 群 算法 在 求解 平台 聚 类 模型 时 具有 优越 性 。 
Rho =0.8 、 信 息 素 总 量 Q=100 、 内 部 负载 权重 a=0.5 、 外 部 


负载 权重 训 =0.5 、 信 息 启发 式 因子 w=08、 期 记 值 启发 式 因子 4。 结束 请 
=3 。 采 用 层次 聚 类 法 和 自 适应 蚁 群 算法 的 输出 平台 聚 类 结 本 文 主 要 解决 了 C2 组 织 中 的 决策 实体 配置 问题 ， 亦 平台 
果 如 表 5、 表 6 所 示 : 聚 类 问题 。 首 先 以 面向 作战 任务 的 形式 从 外 部 和 内 部 两 方面 定 
表 5 传统 层次 聚 类 法 解决 平台 聚 类 问题 的 输出 结果 义 了 决策 实体 的 工作 负载 , 并 以 工作 负载 的 RMS 值 为 目标 函数 
决策 实体 平台 DM 工作 负载 ”工作 负载 RMS 值 建立 了 决策 实体 配置 问题 模型 ， 采 用 基于 全 局 更 新 规则 的 蚁 群 
DMI P,Ps, Po, Pia Py 240 算法 对 该 问题 进行 求解 ， 最 后 通过 仿真 算 例 验证 了 算法 的 适用 
DM; P, ,Po Ps,. Pi 280 性 与 优越 性 。 实 验 结果 表明 ， 合 理 的 算法 参数 设置 能 够 得 到 性 
DM Ps, Py, Pis 200 258.1859 能 较 好 的 聚 类 结果 ， 并 且 该 方法 与 传统 层次 聚 类 法 能 够 得 到 更 
DhM Pa, Pio Pis Pi9 330 小 的 目标 函数 值 ， 且 从 收敛 速度 来 看 该 方法 在 解决 平台 聚 类 问 
DMs Pi, Po, Pia, P17 220 题 时 优 于 传统 蚁 群 算法 和 自 适 应 蚁 群 算法 ， 能 够 较 好 的 解决 决 
表 6 自 适 应 蚁 群 算法 解决 平台 聚 类 问题 的 输出 结果 策 实体 配置 问题 。 在 下 一 步 的 工作 中 ， 将 针对 不 确定 环境 下 的 
决策 实体 平台 DM 工作 负载 工作 负载 RMS 值 决策 实体 配置 方案 优化 调整 问题 进行 研究 。 
DM Ps, P, Pis, Ps. 200 
DM; P, ,Ps, Ps, Po Pyo 260 有 
DM;3 Pio, Pu, Pa P17, Pis, P19 360 257.8759 [1] Bui H, Han X, Mandal S$, et al. Optimization-based decision Support 
DMia Pa, P14 130 algorithm for a team-in-the-loop planning experiment [Cl]// Proc of IEEE 
DMs Ps, Pe, Ps, Pis 280 international Conference on Systems, Man, and Cybernetics, 2009: 
表 4~6 对 比 结果 可 知 ， 采 用 全 局 更 新 规则 的 蚁 群 算法 得 4685-4689. 
到 的 最 终 聚 类 目标 函数 值 为 RMS=256.2421， 层 次 聚 类 法 聚 类 [2] 阳 东 升 , 张 维 明 , 刘 忠 ， 等. 指控 组 织 设计 方法 [M]. 北京 : 国防 工业 
结果 为 RMS=258.1859 ， 自 适应 蚁 群 算 法 聚 类 结果 为 出 版 社 , 2010: 161-189. 
RMS=257.8759。 就 聚 类 结果 而 言 , 采用 全 局 更 新 规则 的 蚁 群 算 。 [3] Han X, Mandal S, Pattipati K R, et al. An optimization-based distributed 
法 得 到 的 聚 类 方案 RMS 值 最 小 ， 自 适应 蚁 群 算法 次 之 ， 传 统 planning algorithm: A blackboard-based collaborative Framework [J]. 


层次 聚 类 法 最 大 。 因 此 蚁 群 算 法 在 解决 平台 聚 类 问题 上 优 于 传 IEEE Trans on Systems, Man, and Cybernetics-Part A: Systems and 
统 层次 聚 类 法 。 而 自 适 应 蚁 群 算法 在 求解 该 问题 时 的 效果 不 如 Humans, 2014, 44 (6): 673-684. 
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